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We present a study of hadronic transitions between T(mS) (m = 4,3, 2) and T{nS) (n — 2, 1) 
resonances based on 347.5 fb~^ of data taken with the BABAR detector at the PEP-II storage 
rings. We report the first observation of T(4S') — > riT{lS) decay with a branching fraction 
B{r{4:S) -^ ?7T(15')) = (1.96 ± OMstat ± O.OQsBsi) x 10"^ and measure the ratio of partial 

widths r{T{4S) -» ■nr{is))/r{r{4S) -> 7r+7r-r(is')) = 2.41 ± 0.40,^* ± 0.12,^,4. We set 90% 

CL upper hmits on the ratios r{r{2S) -^ r^T{lS))/r{T{2S) -^ 7r+7r"r(lS')) < 5.2 x 10"^ and 
r(r(3S') -^ ■nT{lS))/r{r{SS) -^ -k+tv-T{IS)) < 1.9 x lO^^ We also present new measurements 



of the ratios r(r(45') -> 7r+7r-r(2S))/r(r(4S) -^ 7r+7r-r(15')) = 1.16 ± O.lQstat ± O.Usyst 

r(r(3S) -^ TT+Tr-r{2S))/r{r{3S) -^ n+-K-r{is)) = 0.577 ± omGstat ± omOsyst. 



and 



PACS numbers: 14.40.Gx,13.25.Gv 



I. INTRODUCTION 

Hadronic transitions between bound states of heavy 
quarkonia [1] are generally studied using the QCD mul- 
tipole expansion model (QCDME) [2]. This succeeds in 
explaining the relative rates of the ipi'^S) —^ riJ/tp and 
'0(25') —f mrJ/ip transitions and the tttt invariant mass 
distributions in ^(25) ^ ttttJ/V', r{2S) -> 7r7rr(15'), 
T{3S) -^ Tnrr{2S) and the recently observed T(4S') -> 
7r~^Tr~T{lS) decays [3, 4]. Until recently the only feature 
that QCDME could not explain was the dipion invari- 
ant mass distribution in the Y{3S) -^ mTT{lS) tran- 
sition [5], for which a number of possible explanations 
have been proposed [6]. The dipion invariant mass dis- 
tribution in T(4S') — > 7r+7r~T(2S') [3] is also in disagree- 
ment with the QCDME prediction and was not predicted 
either by the alternative explanations proposed for the 
T{3S) -^ 7T^Tr~T{lS). This implies that additional ex- 
perimental input is needed to understand hadronic tran- 
sitions. In QCDME the gluon radiation from a heavy qq 
bound state is calculated in terms of chromo-electric and 
chromo-magnetic fields, in analogy to electromagnctism. 
Transitions between colorless hadrons require the emis- 
sion of at least two gluons. The T{mS) -^ 'imT{nS) 



transitions {m'^Si 



TTTT n Si in spectroscopic nota- 



tion [7]) are ElEl, i.e. transitions where both gluons 
are in an El state. The decays T{mS) -^ ■qT{nS) 
{m^Si -^ rj n^Si) proceed either via E1M2 or MlMl 
transitions; the E1M2 transition is expected to domi- 
nate. The bb system offers unique opportunities: there 
are five known rrv^Si — > tttt n^Si transitions and also 
four kinematically allowed transitions involving an 77 me- 
son. Of the latter only the T{2S) -^ rjT(lS) has been 
recently observed by CLEG [8] , with a branching fraction 
B{T{2S) -> vT{lS)) = {2.ltol ± 0.5) x IQ-*. 

In this paper we present improved measurements of 
the T(4:S) -^ T(nS) transitions, a search for T{mS) -^ 
riT{lS) and new measurements of Y(3S) —^ TT~^Tr~Y(nS) 
and ^(25") ^ 7r+7r-T(lS') partial widths. We also 
measure the ratios of partial widths r(T(?7i5') -^ 
7]r{lS))/T{r{mS) -^ 7r+7r-T(lS')) and T{r{mS) -^ 
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Tr+ IT- r{2S))/T{r{mS) -^ 7r+7r-T(lS')) (m == 3,4), for 
which a number of systematic uncertainties cancel. 

The T{mS) -^ Tr+TT-T{nS) and T{mS) ^ rjTinS) 
transitions, denoted by mS -^ tttt nS and mS -^ 77 nS, 
respectively, are studied by reconstructing the Y{nS) 



mesons via their leptonic decay to fJ.'^ fJ- or e 



+^- 



The 



rj meson is reconstructed via its Tr+Tr^Tr" decay. With 
the choice of this particular 77 decay mode all final states 
contain the same charged particles, resulting in larger 
cancellations of the systematic uncertainties for the ra- 
tios of partial widths. Events where the 77 decays to 77 
are not considered in this work because the €£77 final 
state has a smaller signal to background ratio than the 
i?+£^7r+7r^7r° final state. 



II. DATA SAMPLES AND DETECTOR 

We search for T(4S') hadronic transitions using a sam- 
ple of (383.2 ± 4.2) X 10^ T{4S) decays corresponding 
to an integrated luminosity, L]^, of 347.5 fb~^ acquired 
near the peak of the T(45) resonance ("on-peak", nomi- 
nal center-of-mass energy, ^/s of about 10.58 GeV) with 
the BaBAR detector at the PEP-II asymmetric-energy 
e+e^ storage rings at SLAC. In addition, a data sam- 
ple corresponding to L^Jft = 36.6 fb^^, collected approx- 
imately 40 MeV below the resonance ( "off-peak" ) is used 
to study some of the backgrounds. Decays of T{3S) and 
T{2S) are studied in events recorded "on-peak" and se- 
lected with an initial state radiation (ISR) photon. The 
ISR photon, preferentially emitted at small angle along 
the beam direction, is not required to be detected. 

The BaBAR detector is described in detail elsewhere [9] . 
Charged-particle momenta are measured in a tracking 
system consisting of a five-layer double-sided silicon ver- 
tex tracker (SVT) and a 40-layer central drift chamber 
(DCH), both embedded in a 1.5-T axial magnetic field. 
Charged-particle identification is based on the specific 
energy loss measured in the SVT and DCH, and on a 
measurement of the photons produced in the fused-silica 
bars of the ring- imaging Cherenkov detector (DIRC). A 
CsI(Tl) electromagnetic calorimeter (EMC) is used to de- 
tect and identify photons and electrons, while muons are 
identified in the instrumented flux return of the magnet 
(IFR). 

Simulated Monte Carlo (MC) events are generated us- 
ing the EvtCen package [10]. The angular distribution 
of generated dilepton decays incorporates the T(nS) po- 
larization, while dipion transitions are generated accord- 
ing to phase space. In the simulation of mS -^ rj IS 
we use the angular distribution dictated by the quantum 
numbers for a vector decay to a pseudoscalar and a vec- 
tor. Secondary photon emission is taken into account 
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FIG. 1: Mu vs AM distributions of candidates after the preliminary selection for the fi^ (left) and ee (right) samples. 
Dashed lines delimit the signal boxes for mS -^ tttt nS transitions. The cluster of events in the lower left corner is due to 
r(3S) -> 7r+7r-r(2S') where the r{2S) subsequently decays to T{IS) X. 



in the simulation of T{mS) produced in ISR. Simulated 
events are passed through a detector simulation based on 
GEANT4 [11], and analyzed in the same manner as data. 



III. EVENT SELECTION 

The events of interest have a lepton pair from the de- 
cay of the Y{nS) resonance of invariant mass, Mu, com- 
patible with the known mass values of the T{nS) [12], 
M{nS), and a pair of oppositely charged pions. 

The signature for mS — > tttt nS transition events is an 
invariant mass difference AM — Mt^t^u — Mu compatible 
with the difference of the masses of the two T resonances, 
M{mS) - M{nS), where Mt^ttU is the TT+n^l+i" invari- 
ant mass. 

The mS —^ rj nS events have two additional photons 
from the tt" decay, a 7r+7r~7r° invariant mass, majr, com- 
patible with the known -q mass, M{r]), and an invariant 
mass difference, AM^ = M^^^u — Mu — "13^ compati- 
ble with M{mS) - M{nS) - M{r]), where Ms^u is the 
7r+7r~7r'^£"'"£~ invariant mass. 

The r.m.s. widths of the reconstructed Mu, ra^.^, AM, 
and AM^ distributions are of the order of 75 MeV/c^, 
12 MeV/c^ 7 MeV/c2 and 10 MeV/c^, respectively. 

Events in the data sample with M« within 350 MeV/c'^ 
of the known M{nS) values and AM within 60 MeV/c^ 
of the values expected for any of the mS —^ tttt nS transi- 
tions were not examined until the event selection criteria 
were finalized. Events outside these regions were used to 
understand the background. Simulated MC events were 
used to model the signal. 

Candidate events have at least 4 charged tracks with a 
polar angle 6 within the fiducial volume of the track- 
ing system (0.41< <2.54 rad). Each lepton candi- 
date is required to have a center-of-mass momentum be- 



tween 4.20GeV/c and 5.25GeV/c. At least one of the 
muons of T{nS) — *■ ji'^ IjL~ candidates must be compat- 
ible with the muon hypothesis based on the energy de- 
posited in the EMC and the hit pattern in the IFR along 
the track trajectory. Similarly at least one of the elec- 
trons of T{nS) -^ e+e" candidates must be compati- 
ble with the electron hypothesis based on the energy de- 
posit in the EMC, the ratio of energy in the EMC to 
the track momentum, and the energy loss in the detec- 
tor material. We require Mfj_fj_[Mee] to be within ±200 
[-350,-^200] MeV/c2 of the nominal r{lS) or T{2S) 
mass. The asymmetric cut in the e^e~ sample is due 
to bremsstrahlung, which causes a long tail in the re- 
constructed Mee distribution at low invariant masses and 
that is partially recovered by an algorithm that combines 
the energy of electron tracks with the energy of nearby 
photons. 

Pairs of oppositely charged tracks, not identified as 
electrons and whose Cherenkov angle in the DIRC, when 
measured, is within 3 cr of the value expected for a pion, 
are selected to form a dipion candidate. The dilepton 
and the dipion are constrained to a common vertex and 
the vertex fit is required to have a x^ probability larger 
than 10-3. 

A large fraction of the remaining background is due 
to e+e-7 and fi'^ fi~j events where a photon converts in 
the detector material and the leptons are reconstructed 
as pions. To reduce this background we reject events 
where the opening angle of the charged pion candidates 
in the laboratory reference frame has cos0^+^- > 0.95, 
or where the invariant mass of the charged tracks asso- 
ciated with the pion candidates, calculated assuming the 
e^ mass hypothesis, satisfies rriconv < 50 MeV/c^. The 
distribution of Mu vs AMfor candidate events after the 
preliminary selection is shown in Fig. 1. 

In the case of T{AS) -^ T{nS) transitions the back- 



TABLE I: Selection efficiencies for all studied transitions, sep- 
arately for T{nS) — > fi^ n" and e+e~ as determined by MC 
simulation. For the mS —^ nn nS transitions we quote both 
the efficiency averaged over phase space, eps, and the effec- 
tive efflciency, £e//, calculated according to Eq. 3. 



IV. SIGNAL YIELDS 



Transition 


Selection efficiency (%) 




fifi 


ee 




£PS £eff 


£PS Seff 


25" ^ TTTT IS 


34.46±0.05 36.62±0.08 


ff.f7±0.03 ff.45±0.f4 


3S -^ TTTT IS 


4f.23±0.05 34.f8±0.20 


24.48±0.05 23.96±0.24 


3S -^ TTTT 2S 


f4.76±0.04 f7.2 ±0.6 


^0 


AS -^ TV-K IS 


4f.53±0.23 44.2 ±f.2 


f8.04±0.f8 f9.7 ±2.4 


45 -^ TTTT 2S 


32.69±0.22 30.2 ±0.8 


6.f7±0.f2 7.9 ±3.4 




e 


£■ 


2S^rilS 


8.25±0.09 


^0 


35" ^ J? f 5 


9.42±0.f0 


3.9f±0.06 


■iS ~*rjlS 


f0.07±0.f0 


3.77±0.06 



ground is larger and the expected signal smaller. For this 
reason we further restrict our selection to events where at 
least one of the two charged pions has a transverse mo- 
mentum greater than 100 MeV/c, rriconv > 100 MeV/c^, 
and the polar angle in the laboratory system of the e^ 
from Y(nS) —> e^e~ is larger than 0.7 rad, to reject ra- 
diative Bhabha events. 

Events with at least two candidate photons of E^ > 



50 MeV and invariant mass 110 < 






< 150 MeV/c^ 



are considered to be rj nS candidates if the tt+tt^tt ' in- 
variant mass is within 35 MeV/c^ of the known rj mass. 
To suppress possible cross-feed from the high statistics 
mS —f TTTT nS transitions we require that mS -^ rj nS 
candidates have AM more than 20 MeV/c^ (« 3a) from 
any of the known M{mS) — M{nS) values. 

We select T{2S) and T(3S') states produced via 
ISR, requiring that the momentum of the recon- 
structed £+£^7r+7r^[7r'^] in the center of mass rest frame, 
Pcand^ is within ±150 MeV/c of the expected value of 
(s-M2(mS'))/(2V^). For T(4S') decays p*„„^ is re- 
quired to be < 200 MeV/c. 

The average efficiency for each of the transitions is 
given in Table I. The efficiency for the tt^ tt^ [Tr^]e'^ e^ 
final state is in all cases smaller than for the 
7r"'"7r~[7r°]/i+/i~ final state due to a trigger-level ineffi- 
ciency introduced by the prescaling of Bhabha scattering 
events, whose signature is given by two electrons of large 
invariant mass and no additional charged track of trans- 
verse momentum greater than 250 MeV/c. Because of 
the limited phase-space available in the 2S -^ rj IS and 
35' —f TTTT 2S decays, the momentum of the charged pions 
is always below the threshold, thus the efficiency for these 
two transitions is nearly zero when the Y(nS) decays to 



A. r{mS) 



-r{nS) 



The AM distributions of events in the final sample for 
the mS —f TTTT nS transitions arc shown in Fig. 2. 

We determine the efficiency corrected signal yield for 
the mS -^ TTTT nS transitions without any assumption on 
the angular distributions of the decays. We divide the 
25 -^ TTTT IS and 35 -^ nn 15 samples into 10x6 bins of 
TOtttt and cosOh, where 712,^^^ is the tt+tt" invariant mass 
and 9h is the helicity angle of the 7r+ , defined as the angle 
between the tt"*" direction in the tttt rest frame and the 
TTTT direction in the candidate T{mS) rest frame. The 
35 — > TTTT 25 and 45 -^ tttt nS samples are divided into 
6x4 bins of rriT^^ and cos 0^ . 

The signal yield in each bin is determined by a fit to 
the AM distribution, by maximizing the unbinned ex- 
tended likelihood to the sum of a background probability 
density function (PDF) and a signal PDF. The signal 
PDF is parametrized by a Voigtian function (convolu- 
tion of a Lorentzian with a Gaussian function), that is 
found to describe well the measured AM distribution for 
simulated events. The background is parametrized by a 
linear function. The resolution parameters for the signal 
PDF are fixed to the values determined by the simula- 
tion, thus the free parameters in the fits for bin i are: 
AM]j , the peak position of the signal distribution, Nl^ 
and N^f, , the number of signal and background events, 
and the background shape parameters. The efficiency 
corrected signal yield for each mS —> tttt nS transition is 
then obtained as 



nbins 



Nr.. 



= E ^^./^« 



(1) 



where nbins is the number of bins (60 or 24) and Si is the 
efficiency in each bin determined from MC simulation. 



B. r{mS) -* J7r(lS') 

Figure 3 shows the 7713^ vs AM^ distributions for 
events selected as mS ^ rj IS candidates. The widths of 
the signal boxes have been chosen as « ±3(7 in both vari- 
ables based on MC simulation: Itostt — rn-^| < 35 MeV/c^ 
and |M(m5)-M(15)-M(77)-AM^| < 30 MeV/c^. For 
the 25 ^ 77 15 transition we require AM^ < 30 MeV/c^ 
because the signal for this transition is expected close to 
the kinematic limit. 

The numbers of candidates in the 25 ^ r; 15 and 35 -^ 
T] 15 signal boxes, shown in Table II, are compatible with 
the backgrounds extrapolated from the sidebands defined 
in Fig. 3. Thus, we have no signal for the 25 -^ rj IS 
and 35 ^ 77 15 transitions. 

We observe 56 candidates for the 45 ^ 77 15 transition 
in the "on-peak" data sample, and no candidates in the 
"off-peak" data sample. We test the hypothesis that the 
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FIG. 2: AM distributions of events in the final sample for the mS -^ tttt nS transitions. Data are shown as crosses. The solid 
hnes are the best fit to the data and are only for illustration purposes: they are performed using the signal PDF described in the 
text with resolution parameters fixed to the values determined on MC events. Dashed lines show the background contribution. 



e+e^ ^ rj IS cross section is the same in the "on-peak" 
and "ofF-peak" samples by calculating the binomial prob- 
ability, V, of observing respectively 56 and events for 
a binomial coefficient oi p = L™7(L™* + L^}*^) = 0.905, 
based on the integrated luminosities of the two samples. 
We obtain V — 4: x 10~^ and thus we attribute the ob- 
served r]T{lS) events to T{4:S) decays. 

The event yields for the 45 — > 77 IS* transition in the ee 
and fifx final states are determined by unbinned extended 
maximum likelihood fits to the AAf„ distribution of the 



sample of events in Fig. 3 having m^-^ within 35 MeV/c^ 
of the known r/ mass. The signal PDF is parametrized by 
a Voigtian function, with resolution parameters fixed to 
the values determined from MC events, while the back- 
ground is assumed to be constant. The free parameters 
in the fits are: AM^ sig, the peak position of the sig- 
nal distribution, N^ig and N^kg, the number of signal 
and background events. The efficiency and acceptance 
are determined from MC samples. The fits are shown in 
Fig. 4. The significance, estimated from the likelihood 



TABLE 11: Results for the products of partial widths and branching fractions for the T{mS) hadronic transitions. Ncand is the 
number of candidates in the signal box, Nbck is the number of background events from the fit or estimated from data sidebands 
as described in the text, Ncorr is the efficiency-corrected number of signal events. The first error is statistical, the second is 
systematic. All upper limits are 90%CL. 



Transition 








Our Measurement 


^^ cand 


Nt,,u 


-^^ corr 


Tee (25) X 


B{T{2S) - 


-^ n+n-r{is)) X B{r{is) -^ fl+^i-) 


;meV) 


2582±28±94 


9036 


156±11 


24319±268 


Tee (25') X 


B(r(2S) - 


-> 7r+7r-r(lS')) X S(r(fS) ^ e+e-) 


meV) 


26f8±60±97 


3139 


230±9 


25202±574 


ree(2S') X 


B{r{2S) - 


-> VT{1S)) X B{T{1S) -^ ^+/i') X B{ri -^ 7r+7r-7r°) 


meV) 


<3.1 





2.5±1.1 


< 28 


ree(3S) X 


B{r{3S) - 


-^ 7r+7r-r(is)) X B{r{is) -^ |l+^i-) 


;meV) 


457±8±18 


4198 


207±10 


9945±174 


r^eiSS) X 


B{r{3S) - 


^ TT+TT-r{lS)) X B{r{lS) -^ e+e-) 


;meV) 


441±12±18 


3604 


1234±20 


9821±261 


ree(3S') X 


B{r{3S) - 


-> n+n-r{2S)) X B{r{2S) -^ e+e-) 


;meV) 


206±11±12 


975 


180±21 


4477±241 


ree(3S) X 


B{T{3S) - 


-^ 7ir{lS)) X B{T{1S) -^ ^i+fi-) X B{ri -^ n+n-n°) 


;meV) 


<2.0 


1 


0.8±0.4 


< 41 


ree(3S) X 


B{T{3S) - 


^ vTilS)) X B{T{1S) -> e+e-) x B{r] -^ tt+tt-ttO) 


;meV) 


<9.6 


4 


2.8±0.8 


< 210 


B{r{4:S) - 


-* Tv+n-r{is)) X B{r{is) ^ ^i+^l-) 


;xio-«) 


1.99±0.16±0.07 


687 


378±11 


739± 60 


B(T(4S) - 


-^ 7r+7r-r(lS)) X B{r{lS) -^ e+e-) 


;xio-'') 


1.76±1.05±0.06 


1057 


934±17 


676±397 


B{r{4S) - 


-^ Tv+7v-r{2S)) X B{T{2S) -^ M+p") 


;xio-«) 


1.65±0.21 ±0.11 


377 


204±8 


615± 78 


B{T{4S) - 


-^ 7r+7r-r(2S)) X B(r{2S) -^ e+e") 


;xio-'') 


1.76±1.03 ±0.11 


251 


206±8 


669±392 


B{T{4:S) - 


-vTilS)) 


X BiT{lS) -^ ^+M") X B{ri -^ 7r+7r-7r°) 


;xio-'') 


1.08±0.17±0.05 


40 


0.2±0.4 


387±60 


B{r{4S) - 


-^VT{IS)) 


X B{r{lS) -^ e+e-) X S(r; -> tt+tt-tt") 


;xio-«) 


1.15±0.29±0.05 


16 


0.7±0.6 


424±106 
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FIG. 3: Distributions of m-i-n vs AM,, for the mS ^ rj IS 
transitions studied. Crosses are for the T{1S) -^ e+e- sam- 
ple and dots are for the T{1S) — > l-i+ fJ-- sample. Solid lines 
delimit the signal box region. Dashed lines delimit the side- 
band regions used for background extrapolation. The signal 
box for the 25 ^ r; IS transition (top left) is at the bound- 
ary of the kinematically allowed region of AM,, and only one 
sideband can be defined. 
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FIG. 4: Fits to the AM,, distribution for 4S ^ ri IS candi- 
dates with T{IS) -^ n+ii- (left) and r{lS) -» e+e" (right). 
Data are shown as crosses. The solid lines show the best fit 
to the data. Dashed lines show the background contribution. 



hypothesis, is 11 a and 6.2 a, respectively in the fifi and 
the ee samples. 

The 90% CL upper limits on the signal yields for the 
35" — > ?7 IS* and 25 ^ ?7 IS* transitions are conserva- 
tively estimated from the numbers of events in the signal 
boxes, taking into account the uncertainties in the efh- 
ciencies [13]. The background level in the fi+n^ sample is 
negligible, and background subtraction in the e+e" sam- 
ple, which also has a lower efhciency, would not affect the 
result. 



SYSTEMATIC UNCERTAINTIES 



ratio na ~ -^2 log [C{Nsig)/C{0)] between a fit that in- 
cludes a signal function and a fit with only a background 



We have considered a number of possible sources of 
systematic uncertainties, in addition to the number of 
T{4S) [15] and the calculated luminosity for ISR events. 



TABLE III: Sources of systematic uncertainties on partial widths or branching fractions and ratios of partial widths, separated 
into errors that cancel in ratios, errors due to lepton identification (ID) and invariant mass that are common to all transitions, 
but difi'er for electrons and muons, and errors that are specific to individual decay modes. All errors arc relative and given in 
percent. We also list the corrections applied to account for differences between data and simulation. 



Source 



data/MC 
corr. 



r{2S) -^ 
7r7rr(ig) yrjis) 



r(35) -^ 

ttttTJIS) TnTT{2S) yTjlS) 



r(4S) -> 

7r7rr(lg) 7r7rr(2g) yrjlS) 



Common systematic errors (cancel in all ratios) (%) 



Number of r(4S') 




- 


- 


1.1 


ISR luminosity 




3.0 


3.0 


- 


Tracking 




1.0 


1.0 


1.0 


Selection 




0.3 


0.3 


0.3 


P*cand cut 




0.3 


0.3 


0.3 





Systematic errors associated to lepton identification or invariant 


mass (%) 




Muon ID 


1.025 


0.6 


0.6 






0.6 


M{fi'^fi^) cut 


1.006 


0.2 


0.2 






0.2 


Electron ID 


1.011 


0.7 


0.7 






0.7 


M(e+e") cut 


0.998 


0.5 


0.5 






0.5 



Systematic errors specific to each mode (%) 



7r° efficiency 


1.033 


- 


3.6 


- 


- 


3.6 


- 


- 


3.6 


Acceptance 




0.3 


- 


1.7 


4.7 


- 


2.6 


6.0 


- 


Fitting 




1.6 


1.6 


1.6 


1.6 


1.6 


1.6 


1.6 


1.6 


Total e+e- (%) 




3.7 


- 


4.1 


- 


5.1 


3.5 


6.5 


4.4 


Total n+^i- (%) 




3.7 


5.1 


4.0 


5.9 


5.1 


3.5 


6.5 


4.3 


Total on ratios (%) 




- 


4.3 


- 


5.5 


4.6 


- 


6.9 


5.0 



The uncertainties in charged track and tt^ reconstruction 
efficiencies are determined by a comparison of data and 
MC events on independent control samples. The system- 
atic uncertainties associated with the event selection, the 
cut on p*„„^, the Mu invariant mass cut, and the lepton 
identification criteria are estimated by comparing the effi- 
ciencies determined from MC samples to the correspond- 
ing efficiencies measured with the ISR mS — > tttt nS 
samples in the modes where there are sufficiently high 
statistics and low background to allow the comparison. 
The efficiencies are determined from the numbers of sig- 
nal events which pass or fail any given cut, after all other 
cuts are applied. 

The systematic uncertainties due to the choice of sig- 
nal and background parametcrizations are estimated by 
using different functions or different parameters, and by 
varying the AM or AM^ fit ranges. The uncertainty in 
the acceptance correction for the mS —^ tttt nS tran- 
sitions is determined by the change in the signal yields 
when using different TOtttt and cos 9h binnings. 

The systematic uncertainties from all these sources are 
summarized in Table III for each transition. The total 
systematic uncertainty is estimated by adding in quadra- 
ture all different contributions. We apply correction fac- 
tors to the efficiency determined from MC events, ac- 
counting for differences between data and MC samples 
in the tt*' reconstruction, in lepton identification, and in 
the Mu cut. 



VI. RESULTS 

The products of branching fractions and partial widths 
for each transition are given Table II. They are deter- 
mined from the efficiency-corrected yield in each mode, 
after correcting for small differences between data and 
MC samples and taking into account the number of 
Y{4S) or the equivalent ISR luminosity, /C. For a narrow 
vector resonance produced in ISR 



IC^Ll 



127r^ 



M{mS) 



-Wis,! 



M^{mS) 



(2) 



where the QED "radiator" function W{s^x) is calculated 
to second order following [16-18]. 

Averaging the results from the e+e^ and the /i^/i^ fi- 
nal states, taking into account the common systematic 
errors, and using the world average values of B{rj —^ 
tt+tt-ttO) and B{r{nS) ^ £+£-) [12] we obtain the par- 
tial widths and ratios of partial widths listed in Table IV. 
In this Table, we also compare our results to the values 
expected for each quantity based on previous measure- 
ments of T{mS) widths and branching fractions. The 
measured values of the Y{2S) and T(35) total widths are 
used to derive the theoretical expectations for branching 
fractions from the predicted partial widths in [2]. 

The values of S(r(45) ^ 7r+7r-r(15)) x S(T(1S') ^ 
/i+^-) and B{r{AS) -> 7r+7r-r(2S')) x B{T{2S) -^ 
fi'^IJ,^) supersede our previously reported values based 



10 



TABLE IV: Our measurements for the products and ratios of partial widths and branching fractions of T{mS) hadronic 
transitions, with comparisons to previous measurements and theoretical expectations. We also report the values of the branching 
fractions that are derived from our measurements using world average values for Tee{nS). All upper limits are 90% CL. The 
values of the last seven branching fractions in this Table (reported below the horizontal line) are not independent from the 
values reported above. The values of B{T{4S) -^ TT^TT~T{nS)) from Ref. [12] and indicated with an asterisk are based on our 
previous measurement [3] performed on a subset of the current sample. As discussed in the text, part of the difference in the 
central values is ascribed to a more accurate estimate of the acceptance. 





This work 


PDG [12] 


Prediction 


ree(2S) X 


B{T{2S) -^ n+Tv-r(lS)) 


(cV) 




105.4±1.0±4.2 


115±5 




r(r(25) - 


-^ 77r(i5))/r(r(25) -^ 7r+7r-r(i5)) 


(xlO^ 


-') 


< 5.2 


< 11 


2.5 [2] 


ree(3S) X 


B{T{3S) ^7v+TT-r{lS)) 


(cV) 




18.46±0.27±0.77 


19.8±1.0 




r(r(3S) ^ 


-> 7r+7r-r(2S))/r(r(3S') -^ n+n-T{lS)) 






0.577±0.026±0.060 


0.63±0.14 


0.3 [2] 


r(r(3S) - 


-^ r]T{lS))/r{ri3S) ^ Tv+n-TilS)) 


(xlO^ 


-') 


< 1.9 


<5 


1.7 [2] 


B{r{4S) - 


-^ 7V + TT-r{lS)) 


(xlO" 


-') 


0.800±0.064±0.027 


0.90±0.15(*' 


- 


r(r(45') - 


^ 7r+7r-r(25))/r(r(4S) ^ 7r+7r-r(lS)) 






1.16±0.16±0.14 




- 


r(r(4S) - 


-^ 77r(i5))/r(r(45) ^ ^+^-r(i5)) 






2.41±0.40±0.12 


- 


- 


B{T{2S) - 


-^ 7r+7r-r(15)) 


(%) 




17.22±0.17±0.75 


18.8±0.6 


27±2 [2] 


B{T{2S) - 


- VT{1S)) 


(xlO- 


-4) 


<9 


<20 


8.1±0.8 [14] 


B{r{3S) - 


-^ 7v+n-r{lS)) 


(%) 




4.17±0.06±0.19 


4.48±0.21 


3.3±0.3 [2] 


B{r{3S) - 


-^ n+'K-r{2S)) 


(%) 




2.40±0.10±0.26 


2.8±0.6 


l.OiO.l [2] 


B{r{3S) - 


- riYilS)) 


(xlO- 


-4) 


< 8 


< 22 


6.7±0.7 [14] 


s(r(4S) - 


^ Tv+n-T{2S)) 


(xlO" 


-') 


0.86±0.11±0.07 


0.88±0.19(*' 


- 


s(r(4S) ^ 


^VT{1S)) 


(xlO- 


-') 


1.96±0.06±0.09 


- 


- 



on a fraction of the current sample [3] . Part of the differ- 
ence in the central values is due to the different methods 
used to determine the acceptance, which was calculated 
in our previous paper assuming a phase-space distribu- 
tion in the T(4S') -^ 'K^iT~T{nS) decay. The efficiency is 
not uniform over the Dalitz plot, thus the impact on the 
central value between the two methods depends on the 
angular distributions peculiar of each transition. The 
difference can be estimated by comparing the value of 
the phase-space averaged efficiencies eps, and the effec- 
tive efficiencies £eff calculated from the observed event 
yields in each region of the Dalitz plot 



^e// 



Enbins T»r 



sig 



Enbins ^i i 
i=l ^^ sig 1^1 



(3) 



Notice that the uncertainty in the calculated effective 
efficiency is due to the statistical uncertainty in the event 
yield. As shown in Table I the effective efficiency for 



T{AS) 



T~^ IT 



T{1S), when the T{1S) decays to /i+/x 



IS ~ 7% larger than the value estimated using a phase- 
space distribution. Accounting for this difference, the 
results presented here are statistically compatible with 
the ones previously reported. 

From our result we derive new values for 
B{T{'iS,2S) -> TT+TT-T{1S)) that are of comparable 
precision to the previous world averages, and compatible 
with them. The value of B{T{5S) ^ 7r+7r-r(2S')) 
derived from our measurement has an error that is 
smaller than the current world average. 



VII. CONCLUSIONS 



We have presented a study of hadronic transitions be- 
tween the T states: new measurements of the branch- 
ing fractions S(T(4S') ^ 7r+7r-T(lS', 2S')), S(T(3S') ^ 
7r+7r^T(2S')) which have smaller errors than current 
world averages, and new measurements oiB{T{3S, 2S) -^ 
TT^TT^T{1S)) whose precision is comparable to present 
world averages. We have also presented measure- 
ments of the ratios of partial widths T{T{mS) -^ 
TT+Tr-r{2S))/T{r{mS) ^ 7r+7r-T(15)) (m = 3,4) 
where a number of systematic uncertainties cancel. Our 
results for the branching fractions of the T{2S) and 
Y{3S) -^ riT{lS) transitions represent improvements 
over the current published upper limits, and are compat- 
ible with the recent results from CLEO [8]: B{T{2S) -> 
7?T(15)) == (2.ll:g;^±0.5) x IQ-^, BiT{3S) ^ 7?T(15)) < 
2.9 X lO"'^ at 90% CL. 

We observe a significant number of r]Y(lS) candi- 
dates at the formation energy of the T(4S'). We can 
exclude the hypothesis that they are due to continuum 
e+e" -^ riT{lS) with a probability of 99.6% and we 
attribute them to T{AS) decays. The branching frac- 
tion for the T{4:S) -^ riT{lS) decay is larger than the 
branching fraction for T(4S') -^ Tr~^TT~T(lS), which is 
unexpected when compared to all other known charmo- 
nium and bottomonium transitions. There are no pre- 
dictions for this specific decay mode. In the QCDME 
calculation for hadronic transitions, the effect of the 
nodes in the wave functions in the overlap integrals be- 
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tween the initial and final states and the intermediate 
states can be large for radial excitations. But even that 
should not significantly affect the ratio of partial widths 

r(r(4S') ^ 'nr{isj)/r{r{AS) -> TT+n-r{is)), at least 

if the T(4S') ^ ?7T(15) transition is E1M2 [2]. It is 
possible that accidental cancellations suppress the E1M2 
term with respect to MlMl, or perhaps QCDME be- 
comes unreliable for higher gluon momenta. These re- 
sults, together with the recent CLEO measurement of 
the matrix elements in T{3S) — > 7r+7r~T(lS', 25') and 
T{2S) -^ 7r+7r~T(lS') transitions [19], could provide a 
tool to understand the hadronic transitions better. 
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